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Automated  DNA  Sizing  by  Atomic  Force  Microscopy 
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j---.)  NA  sizing  is  one  of  the  most  widely  used  analytical  methods  in  molecular 
biology  and  biochemistry,  and  it  is  a  core  method  for  genomic  analysis.  Many  research, 
health  care,  and  forensic  applications  rely  on  DNA  sizing.  Given  the  volume  of  DNA 
sizing  done  and  the  associated  cost,  an  increasing  need  exists  for  extremely  rapid,  high- 
throughput,  and  inexpensive  techniques.  Essentially  all  DNA  sizing  today  is  performed 
by  gel  electrophoresis,  although  multiple  efforts  are  under  way  to  develop  more  effective 
technologies  based  on  entirely  new  approaches.  One  such  effort  is  described  in  this 
article:  the  use  of  atomic  force  microscopy  in  combination  with  pattern  recognition 
software  to  perform  DNA  sizing.  Initial  results  show  that  this  approach  can  be 
compared  with  gel  electrophoresis  for  some  applications.  Future  work  aims  at 
developing  a  practical,  high-throughput  DNA  sizing  system.  (Keywords:  Atomic  force 
microscopy,  DNA  sizing.  Image  processing.) 


INTRODUCTION 

Purpose  of  Sking  DNA 

Why  is  DNA  sizing  important?  Combined  with 
other  information,  the  size  of  a  DNA  fragment  can 
contain  information  about  the  sequence  within  the 
fragment,  or  the  DNA  from  which  the  fragment  was 
derived.  This  information  can  be  used  in  a  wide  range 
of  applications,  such  as  constmction  of  physical  ge¬ 
nome  maps  and  genotyping.  Physical  maps  of  DNA 
clones  provide  key  information  for  completing  the 
sequence  of  large  genomes  (e.g.,  the  Human  Genome 
Project).  Physical  map  constmction  is  also  important 
for  positional  cloning,  a  process  of  isolating  genes  that 
have  only  been  identified  genetically.  One  type  of 
physical  map,  called  a  restriction  map,  represents  the 


locations  along  the  DNA  cut  by  restriction  enzymes. 
Genotyping  is  the  process  of  characterizing  genetic 
differences  throughout  the  genome  among  different 
people.  These  differences  are  used  to  follow  inheritance 
of  genes,  to  identify  people  having  variants  of  genes 
(alleles)  that  predispose  to  specific  diseases,  and  to 
clearly  identify  different  people,  such  as  in  forensic 
applications. 

The  genomes  of  most  organisms  are  too  large  to 
examine  as  a  whole,  so  DNA  is  divided  into  manageable¬ 
sized  fragments.  For  example,  sequencing  the  entire 
human  genome  of  approximately  3  billion  base  pairs 
(bp)  requires  that  much  smaller  sequences  (=500  bp) 
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be  determined  at  a  timed  The  process  of  breaking  the 
DNA  into  manageable  pieces  and  sorting  those  pieces 
requires  DNA  sizing. 

The  information  produced  from  DNA  sizing  mainly 
depends  on  how  the  fragments  are  generated,  although 
biochemical  techniques  are  used  for  most  sizing  appli¬ 
cations.  Some  common  methods  are  synthesis  of  com¬ 
plementary  strands  from  a  template  DNA,  cutting 
DNA  at  specific  sites  by  restriction  enzymes,  or  repli¬ 
cating  a  small  portion  of  the  DNA  using  a  process 
called  PCR  (polymerase  chain  reaction). 

Approaches  to  DNA  Sizing 

Today,  gel  electrophoresis  methods  are  used  fre¬ 
quently  for  DNA  sizing.  By  combining  the  gel  type  and 
size  with  the  applied  electric  field  strength,  gel  electro¬ 
phoresis  can  be  used  to  size  DNA  molecules  from  a  few 
bases  to  millions  of  bases.^  This  method  is  relatively 
simple,  the  equipment  needed  is  inexpensive,  and  the 
results  are  widely  understood  and  accepted.  Hence,  for 
laboratory  research,  the  approach  is  very  useful.  Gel 
electrophoresis  has  become  extremely  sophisticated  in 
recent  years,^'^  but  in  principle  the  current  technology 
is  the  result  of  incremental  improvements  to  the  basic 
approach  that  is  over  40  years  old.  Its  primary  limita¬ 
tions  are  the  slow  speed  (running  one  gel  often  takes 
2  h  or  more)  and  the  need  for  large  sample  amounts  of 
DNA  (typically  hundreds  of  nanograms) — both  factors 
that  directly  affect  cost.  Thus,  alternative  approaches 
aim  at  increasing  speed  and  reducing  the  amount  of 
sample  required. 

Of  the  alternative  approaches  now  emerging,  the  use 
of  optical  microscopy  for  genomic  analysis  is  the  most 
mature.^^  Optical  mapping  uses  light  microscopy  to 
determine  the  length  of  fluorescently  stained  DNA 
restriction  fragments.  This  method  has  advantages  over 
the  gel-based  method  in  throughput,  resolution,  safety, 
and  cost.  However,  optical  microscopy  is  not  very  ef¬ 
fective  for  smaller  DNA  fragments;  the  resolution  is 
diffraction  limited  at  about  200  nm,  which  is  approx¬ 
imately  600  bp. 

Other  approaches  to  DNA  sizing  include  flow  cy¬ 
tometry,  mass  spectrometry,  electron  microscopy,  elec¬ 
trophysiology,  and  chip-based  technologies.  Flow  cy¬ 
tometry,  which  sorts  DNA  molecules  based  on  labeling 
with  fluorescent  markers,  has  been  used  extensively  to 
size  and  separate  larger  DNA  fragments  such  as  chro¬ 
mosomes.  At  this  time  its  lower  limit  has  been  im¬ 
proved  to  about  1500  bp,®’^  even  larger  than  the  lower 
limit  for  optical  microscopy.  Mass  spectrometry  is  also 
emerging  as  a  tool  for  genotyping  and  possibly  DNA 
sizing.^°’^^  Today,  however,  it  cannot  easily  be  used  on 
DNA  molecules  larger  than  100  bp  in  length.  Chip- 
based  technologies  have  taken  over  certain  applica¬ 
tions  for  which  DNA  sizing  by  gel  electrophoresis  was 


previously  used,  but  do  not  report  fragment  size.  Elec¬ 
tron  microscopy  is  widely  used  for  visualizing  the  struc¬ 
ture  of  small  objects,  including  DNA.  However,  the 
relatively  cumbersome  sample  preparation  and  opera¬ 
tion  of  the  instrument  have  prevented  its  use  in  high- 
throughput  applications.  A  recently  proposed  approach 
for  DNA  sizing  is  the  use  of  ion  channels  and  electro- 
physiology.^^  In  this  approach,  the  occlusion  of  the 
chaimel  by  DNA  as  it  passes  through  results  in  changes 
in  current  flow  through  the  channel.  This  blockage  of 
the  channel  is  quantified  using  electrophysiological 
methods  and  is  related  to  the  length  of  the  DNA. 

The  atomic  force  microscope  obtains  height  images 
by  scanning  a  surface  with  an  extremely  small  tip  at  the 
end  of  a  cantilever.  The  height  is  determined  at  posi¬ 
tions  along  the  scan  using  voltage  measurements  sam¬ 
pled  from  a  split  segment  photodiode  array  aligned  with 
a  laser  beam  deflected  off  the  top  of  the  cantilever,  as 
shown  in  Fig.  1 .  Optimally,  the  atomic  force  microscope 
can  resolve  single  atoms  (0.2-nm  resolution). Under 
more  typical  imaging  conditions  the  resolution  is  2  to 
20  nm,  which  is  still  more  than  an  order  of  magnitude 
better  than  optical  microscopy  (>200  nm).  The  next 
section  describes  the  merits  of  DNA  sizing  using  atomic 
force  microscopy  (AFM). 

ATOMIC  FORCE  MICROSCOPY 

It  is  well  established  that  AFM  can  easily  visualize 
nucleic  acids  adsorbed  to  a  surface,  and  that  these 
images  reflect  the  size  of  the  DNA  molecule.*'’”^’  Figure 
2  shows  a  typical  AFM  image  of  DNA  fragments  on  a 
mica  substrate.  The  imaging  of  DNA  can  be  done  in 
liquids  or  on  air-dried  samples,  making  AFM  very  prac¬ 
tical.  Although  no  fundamental  upper  size  limit  exists 
for  objects  that  can  be  imaged  by  AFM,  issues  such  as 
bandwidth  and  sample  preparation  make  AFM  more 


Figure  1.  An  atomic  force  microscope  scans  a  surface  with  an 
extremeiy  smaii  probe  tip  at  the  end  of  a  cantilever.  Voltage  measure¬ 
ments  from  a  laser  beam  deflected  off  the  top  of  the  cantilever  to  a 
split  segment  photodiode  are  used  to  create  height  images. 
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Figure  2.  An  AFM  image  of  DNA  on  mica. 


suitable  for  DNA  smaller  than  10,000  bp.  We  estimate 
the  practical  lower  limit  to  be  less  than  50  bp,  depending 
on  the  imaging  resolution.  At  smaller  sizes  of  DNA,  the 
convolution  of  the  image  with  the  cantilever  tip  shape 
obscures  the  shape  (and  thus  the  length)  of  the  molecule. 

Hence,  the  question  of  whether  AFM  becomes  a 
useful  DNA  sizing  technology  lies  not  in  its  ability  to  size, 
but  rather  whether  a  practical  and  useful  implementa¬ 
tion  can  be  developed.  This  question  of  implementation 
is  at  the  core  of  our  present  efforts.  AFM  has  been  used 
extensively  for  quality  control  of  semiconductors.  These 
systems  use  atomic  force  microscopes  that  run  24  h  per 
day  with  automated  sample  exchange,  tip  positioning  (to 
better  than  l-)Jm  accuracy),  tip  exchange,  and  data 
validation.  Thus,  many  implementation  issues  for  auto¬ 
mating  such  a  microscope  have  been  solved.  Others  yet 
to  be  solved  are  sample  deposition,  high-resolution  im¬ 
aging  in  a  noisy  environment,  and  automated  pattern 
recognition  for  DNA  size  determination. 

In  general  terms,  AFM  falls  into  the  category  of 
single  molecule  technologies.  One  of  the  strongest 
trends  in  analytical  chemistry  in  the  past  several  de¬ 
cades  is  the  miniaturization  of  technology,  which  is 
driven  by  considerations  such  as  reduced  sample  size, 
reagent  cost,  and  labor  and  analysis  time.  At  the  low 
end  of  the  road  to  miniaturization  are  the  single¬ 
molecule  detection  techniques,  including  AFM.  Thus, 
AFM  capabilities  are  consistent  with  current  trends  in 
analytical  technology. 

An  additional  advantage  of  AFM  is  the  high  signal- 
to-noise  ratio,  which  allows  direct  visualization  of 
individual  DNA  or  protein  molecules  without  contrast 
enhancing  agents.^'*”*®  Further,  AFM  has  the  potential 
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to  directly  detect  details  of  DNA  such  as  supercoiled 
geometries, ^^kinks,^°  and  varied  helical  repeats  of  sin¬ 
gle  molecules, and  to  directly  map  the  specificity  and 
structure  of  DNA^binding  proteins  bound  to  mole- 
cules,^^-^^  thereby  providing  additional  types  of  valu¬ 
able  information.  These  added  capabilities  can  have 
considerable  advantages  over  other  approaches  and  are 
discussed  later  in  this  article. 

General  Approach 

The  present  approach  is  founded  on  two  premises: 
First,  sizing  of  DNA  will  continue  to  be  a  major  ana¬ 
lytical  method  used  in  genomic  analysis  for  the  foresee¬ 
able  future.  Thus,  the  proposed  work  will  not  be  “or¬ 
phaned”  by  shifts  in  specific  applications.  Second, 
AFM  currently  has  a  higher  signal-to-noise  ratio  than 
any  other  analytical  technology  that  can  be  applied  to 
nucleic  acids.  This  has  two  primary  consequences: 

1 .  AFM  can  determine  the  size  of  a  single  piece  of  DNA 
whose  length  is  near  or  below  the  optical  diffraction 
limit. 

2.  AFM  can  directly  detect  small  proteins  and  other 
markers  bound  to  a  single  piece  of  DNA  without 
contrast  enhancement. 

The  general  strategy  proposed  is  shown  in  Fig.  3  as 
a  comparison  to  gel  electrophoresis.  A  DNA  sample  in 
solution  is  adsorbed  to  a  solid  support  and  imaged  by 
AFM;  the  lengths  of  the  DNA  molecules  in  the  sample 
are  determined  by  automated  image  processing  soft¬ 
ware,  resulting  in  a  histogram  of  sizes  essentially  iden¬ 
tical  to  that  produced  by  the  scan  of  an  electrophoresis 
gel.  We  call  this  approach  solid-state  DNA  sizing 
(SSDS).^®  The  results  indicate  that  SSDS  clearly  works 
in  principle  and  on  a  limited  scale  in  practice.  As 
mentioned  earlier,  the  core  issue  is  implementing  this 
approach  and  making  it  practical. 

Automatic  Sizing  Software 

To  determine  rapidly  and  accurately  the  distribution 
of  DNA  lengths  in  AFM  images,  we  have  developed 
specialized  software  written  in  the  Interactive  Data 
Language.  This  language  is  a  high-level  interactive  pro¬ 
gramming  environment  designed  to  handle  and  manip¬ 
ulate  large  arrays  of  data  and  has  many  advanced  pro¬ 
cessing  took.  The  automatic  sizing  software  was 
developed  to  be  platform  independent  and  has  a  user- 
friendly  graphical  interface.  The  drawback  of  this  lan¬ 
guage  is  that  it  is  slower  than  typical  compiled  computer 
languages  such  as  C/C++.  However,  the  processing  time 
per  image  is  still  less  than  60  s,  which  allows  more  than 
1500  images  per  day  to  be  processed  on  an  inexpensive 
personal  computer. 

The  input  to  the  software  is  the  raw  AFM  image  data 
of  DNA  adsorbed  onto  a  substrate.  The  software 
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Figure  3.  A  schematic  for  DNA  sizing  based  on  AFM  imaging 
proposed  here  compared  with  agarose  gel  electrophoresis.  DNA 
molecules  are  adsorbed  onto  a  solid  support  (mica).  The  adsorbed 
DNA  is  imaged  by  AFM,  and  the  fragment  lengths  are  automatically 
determined  using  the  developed  sizing  software.  The  result  is  a 
histogram  of  fragment  lengths  similar  to  the  scan  of  a  DNA  sample 
separated  on  an  agarose  gel.  (Reprinted  from  Ref.  26  by  permis¬ 
sion.  ©1998  American  Chemical  Society.) 


executes  a  series  of  conventional  image  processing  steps 
(Fig.  4).  The  end  result  after  a  pixel-to-pixel  length 
determination  of  each  fragment  is  a  histogram  of  sizes 
similar  to  a  scan  of  an  electrophoresis  gel  of  separated 
DNA  (for  details  of  the  algorithm,  see  Ref.  27).  The 
most  critical  data  parameters  for  this  software  are  that 
the  image  background  is  kept  to  a  reasonable  level  and 
the  density  of  fragments  is  not  so  high  that  frequent 
overlaps  are  present.  Any  fragments  that  overlap  or 
touch  the  edge  of  the  image  are  excluded  by  the  sizing 
software  to  avoid  erroneous  results. 

Experiments  and  Results 

The  accuracy,  precision,  resolution,  and  sensitivity 
of  AFM-based  DNA  sizing  for  small  DNA  fragments 
were  evaluated.  Single  samples  were  prepared,  imaged, 
and  processed  with  the  image  processing  software. 
Detailed  descriptions  of  these  tests  and  results  are 


published  elsewhere.^^^®  A  summary  is  presented  here. 
The  length  of  the  DNA  in  these  tests  is  given  both  in 
base  pairs  and  nanometers,  and  the  relative  scale  for 
these  DNA  is  1  bp  equals  0.338  nm. 

Accuracy  is  determined  by  testing  a  series  of  DNA 
samples  with  known  size.  Typically  the  measured  length 
was  within  15%  of  the  expected  length  based  on  the 
known  size.  The  precision  was  assessed  by  the  width  of 
the  distribution  peak  at  half  height.  This  value  was 
typically  less  than  10%  of  the  fragment  length  for  frag¬ 
ments  larger  than  200  bp.  The  sizes  obtained  by  hand 
tracing  the  DNA  in  the  images  are  similar  to  those 
obtained  with  the  automated  sizing  program  (Fig.  5). 
The  hand  tracing  was  accomplished  by  using  NIH 
Image  software  and  interactively  tracing  over  the  DNA 
fragment.  Results  from  hand  tracing  can  vary  depending 
on  the  way  in  which  the  end  points  of  the  molecule  are 
determined;  the  sizing  software  produces  the  same  re¬ 
sults  each  time. 

The  precision  also  gives  one  measurement  of  the 
resolution.  A  more  practical  assessment  of  resolution  is 
provided  by  examining  known  mixtures  of  DNA  frag¬ 
ments.  Figure  6  shows  the  results  from  an  experiment 
where  a  200-bp  fragment  can  be  distinguished  from  a 
250-bp  fragment.  In  another  experiment  with  a  sample 
containing  fragments  from  100  to  2000  bp,  the  different 
fragments  can  be  easily  separated  in  the  histogram  of 
size  distributions  shown  in  Fig.  7.  The  longer  DNA 
fragments  in  these  tests  are  less  well  represented  both 
in  the  histograms  and  in  the  images  that  were  processed 
because  the  adsorption  kinetics  of  longer  fragments  are 
relatively  slow. 

For  a  practical  test  application,  a  blind  comparison 
with  agarose  gel  electrophoresis  for  screening  a  set  of 
PI  artificial  chromosomes  (PACs)  was  completed.  A 
PAC  library  (BAC  PAC  Resources,  Roswell  Park  Can¬ 
cer  Institute,  Buffalo,  New  York)  was  screened  to  gen¬ 
erate  the  PAC  set  examined.  A  single  yeast  artificial 
chromosome  (YAC)  of  350  kb  was  isolated  from  a 
pulsed  field  gel,  radioactively  labeled,  and  hybridized  to 
PAC  library  filters.  The  resulting  20  positive  clones 
were  screened  for  the  presence  of  known  markers  on  the 
YAC  using  PAC  DNA  as  templates  for  either  15-  or  35- 
cycle  PCR  reactions.  The  same  set  of  PCR  products  was 
examined  by  the  AFM-based  sizing  method.  The  results 
from  SSDS  were  in  excellent  agreement  with  the  gel- 
based  results  (Fig.  8). 

The  error  rate  remains  to  be  formally  established, 
although  we  estimate  the  present  error  rate  to  be  be¬ 
tween  0.1  and  1%.  Errors  typically  occur  when  the 
sample  preparation  fails  in  some  way,  so  that  no  sample 
is  present  on  the  substrates.  The  time  required  for  SSDS 
is  several  minutes  per  sample.  For  a  small  number  of 
samples  this  technique  is  better  than  gel  electrophore¬ 
sis,  which  often  requires  hours  to  perform.  However,  gel 
electrophoresis  allows  for  simultaneous  processing  of 
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Figure  4.  Examples  of  intermediates  from  the  automated  DNA  sizing  software,  (a)  AFM  image,  (b)  After  adaptive  thresholding  based  on 
the  optimum  threshold  obtained  by  Gaussian  fitting  of  both  foreground  and  background  gray  scale  data,  and  followed  by  smoothing  the 
binary  image  with  a  3  x  3  pixel  average,  (c)  After  removing  edge-touching  fragments  and  very  small  artifacts,  (d)  After  pruning  and  thinning 
by  using  a  fast  parallel  thinning  algorithm,  the  pixels  are  removed  only  if  a  set  of  conditions  for  each  pixel’s  neighborhood  is  met,  until  no 
further  pixels  satisfy  the  conditions  for  removal,  (e)  After  removing  overlapping  fragments,  returning  end-pixels  that  were  deleted  by 
thinning,  and  then  calculating  the  length  in  nanometers  of  each  remaining  fragment,  (f)  Histogram  of  the  calculated  lengths. 
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Figure  5.  Comparison  of  automated  sizing  to  hand  tracing  of  the 
same  image  of  a  sample  containing  fragments  250  bp  long, 
showing  the  accuracy  and  consistency. 
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Figure  6.  Empirical  evaluation  of  the  resolution  using  a  sample  of 
an  equimolar  mixture  of  200-  and  250-bp  DNA,  showing  two 
resolved  peaks  near  the  actual  sizes  in  the  histogram. 
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Figure  7.  Length  distribution  of  common  DNA  reference  sizes 
mixed  in  a  singie  sampie  containing  100-,  200-,  400-,  800-,  1200-, 
and  2000-bp  fragments,  indicating  the  good  sizing  capabiiity  over 
a  iarge  range. 


several  dozen  samples.  Thus,  the  time  per  sample  is 
similar  for  the  two  approaches.  (In  one  set  of  experi¬ 
ments,  one  person  completed  manual  processing  of  60 
PCR  samples  in  about  4  h.)  Many  steps  of  this  process 
are  amenable  to  automation  and  should  provide  rapid 
throughput.  The  emergence  of  parallel  tip  arrays  may 
further  improve  the  throughput  in  SSDS. 

Other  Capabilities  of  AFM 

The  high  signal-to-noise  ratio  of  AFM  and  the  re¬ 
sulting  ability  to  detect  proteins  bound  to  the  DNA,  or 
other  sequence-specific  modifications,  lead  to  signifi¬ 
cant  potential  advantages  over  gel  electrophoresis.  In 
particular  it  now  becomes  possible  to  “decorate”  the 
DNA  in  a  way  that  increases  the  specific  information 
obtained  about  the  molecule.  In  fact,  decorating  DNA 
is  analogous  to  bar  coding  because  it  specifies  the  identity 
of  a  piece  of  DNA  with  a  certain  level  of  confidence. 

Decorating  DNA  is  relevant  to  some  sizing  applica¬ 
tions  such  as  restriction  mapping,  but  it  is  more  power¬ 
ful  than  restriction  mapping  as  currently  implemented. 
Restriction  mapping  locates  sequence-specific  positions 


along  a  piece  of  DNA  that  are  cut 
by  restriction  enzymes.  The  proce¬ 
dure  involves  first  sizing  the  DNA 
fragments  resulting  from  a  piece  of 
DNA  cut  by  enzymes.  Then  the 
order  of  the  fragments  along  the 
original  piece  of  DNA  must  be  re¬ 
solved,  which  is  an  inefficient  and 
redundant  process.  The  high  con¬ 
trast  of  AFM  allows  bound  proteins, 
such  as  “stalled”  restriction  en¬ 
zymes,  to  be  visualized  on  the 
DNA,  thereby  allowing  restriction 


maps  to  be  determined  directly  by 
the  distance  between  the  bound 
proteins  on  the  DNA.^’^®’^°  This 
ability  to  directly  determine  the 
sequence-specific  decoration  pat¬ 
tern  eliminates  the  step  required  to 
resolve  the  order  of  the  fragments. 

Figure  8.  Comparison  of  AFM-based  siz¬ 
ing  with  agarose  gel  electrophoresis  for 
screening  370-bp  DNA  products.  Results 
of  thetwo  methods  were  in  excellent  agree¬ 
ment.  (a)  Number  of  DNA  fragments  de¬ 
tected  within  a  30-nm  window  of  the  ex¬ 
pected  size  for  each  of  the  20  samples 
using  AFM-based  sizing,  (b)  Section  of  the 
gel  along  with  a  manually  placed  scan  line 
across  the  DNA  size  to  be  screened  (370 
bp),  (c)  Plot  of  the  intensity  of  the  line  scan 
indicating  the  relative  numberof  DNA  frag¬ 
ments  in  each  of  the  20  samples.  (Re¬ 
printed  from  Ref.  26  by  permission.  ©  1 998 
American  Chemical  Society.) 
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In  addition,  more  types  of  sequence-specific  decorations 
exist  that  can  be  used  with  AFM  than  just  the  restric¬ 
tion  enzymes  used  for  restriction  mapping.  Therefore, 
developing  and  implementing  this  decoration  capabil¬ 
ity  is  also  a  high  priority  for  our  future  efforts. 


CONCLUSION 

The  conditions  for  AFM-based  DNA  sizing  estab¬ 
lished  here  provide  a  basis  for  development  of  a  highly 
automated,  high-throughput  assessment  of  DNA  sam¬ 
ples  with  a  variety  of  applications  in  research,  health 
care,  and  forensics.  The  potential  for  reduced  time, 
smaller  sample  amounts,  and  lower  reagent  costs,  as  well 
as  the  available  technologies  for  automation,  suggest 
that  AFM  represents  a  competitive  approach  for  DNA 
sizing  and  other  types  of  DNA  analyses. 
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